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Subdiffusive transport in tilted washboard potentials is studied within the fractional Fokker- 
Planck equation approach, using the associated continuous time random walk (CTRW) framework. 
The scaled subvelocity is shown to obey a universal law, assuming the form of a stationary Levy- 
stable distribution. The latter is defined by the index of subdiffusion a and the mean subvelocity 
only, but interestingly depends neither on the bias strength nor on the specific form of the potential. 
These scaled, universal subvelocity fluctuations emerge due to the weak ergodicity breaking and 
are vanishing in the limit of normal diffusion. The results of the analytical heuristic theory are 
corroborated by Monte Carlo simulations of the underlying CTRW. 
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A process of directed motion, for example the motion 
of a Brownian particle under influence of a constant force, 
can be characterized by its mean velocity v. The mean 
velocity is measured using a ruler and a stopwatch in one 
of two different setups: One can measure the distances L 
covered within a fixed time interval t, or, like it is done in 
sport competitions, fix the distance L and measure the 
time intervals necessary to cover it. Thus, one can distin- 
guish between the the fixed time (FT) velocities and the 
time-of-flight (TOF) velocities. In any case, in "normal" 
situation the typical time necessary to overcome the dis- 
tance L grows on the average linearly with L in the TOF 
setup, or the typical distance covered during the time t 
grows on the average linearly with t in the FT setup. For 
our Brownian particle moving under the influence of the 
constant force F both setups give the values oi L/t which 
in the limit of i — s- oo for the FT measurement or L — > oo 
for the TOF measurements reach the same sharp value 

V. 

For the case of biased anomalous diffusion (subdiffu- 
sion) the situation drastically differs. In what follows the 
subdiffusive motion x(t) is modeled by a continuous time 
random walk (CTRW) with the waiting time probability 
density (WTD) on sites following (for t ^ t) a. power 
law 



^P{t) ^ C{t/T)- 



(1) 



with a diverging mean waiting time, i.e. with < a < 1. 
In ([1]), r is a characteristic time scale and the prefactor 
c — [TQ;r(l — a)]^^ is introduced for simplicity of further 
calculations. 

For example, charge transport processes in disordered, 
amorphous media can be subdiffusive due to a trap- 
like transport mechanism with a similar to ([T]) trap- 
ping time distribution P, 0, 0, 0j S Q; an approxima- 
tion which can be justified for samples of macroscopic, 
but finite size L Q. The corresponding averaged cur- 
rent J{t) oc d{Sx{t)) /dt, Sx{t) = x{t) — x{0), is a tran- 



sient, decaying to zero quantity [3, [^. However one 
can define an ensemWe-averaged mean subvelocity Va = 
which is a quasi-stationary quan- 
tity for a sufficiently large L (neglecting finite size effects, 
i.e. L ^ OO when assuming limit i — > oo). In case of de- 
caying photocurrent experiments in thin amorphous films 
3], one can define an analogue of subvelocity, namely 
anomalous current as Jq — J{t')dt' /t" = Q{t)/t°', 
where Q{t) is the transferred charge. It will first be quasi- 
stationary and then decay anyway due to the finite size, 
edge effects. 

The absence of a mean trapping time leads to the 
(weak) ergodicity breaking [9, 10] in the relevant trans- 
port processes. In particular, the mean subvelocity of 
individual particles (before ensemble averaging) is a ran- 
dom quantity as the time and ensemble averages are not 
equivalent. This is just like the diffusion coefficient of in- 
dividual particles is a random quantity within our setup 
pll . [T^ . We shall consider a CTRW in a tilted periodic 
potential. In the continuous space limit, it is described 
by the fractional Fokkcr-Planck equation (FFPE) ^, 7,, 8] 



d°'P{x,t) 



d_ 

dx 



e-™|.eW-)p(x,t) 
ox 



(2) 



which we write down here in the form with the fractional 
Caputo derivative a"P(a;,i)/9r = (l/r(l-a)) Jldt'[t- 
t']-°'dP{x,t')/dt' In Eq. @, U{x) = V{x) - Fx, 

where V{x-\-l) = V{x) is a periodic potential with period 
I, and F > is the biasing force; (3 = 1/(A;bT) is the in- 
verse temperature, and the subdiffusion coefficient. 
The self-averaging in time does not occur 0] and mean 
subvelocity remains a random variable even in the strict 
limit t ^ OO. To find the corresponding probability dis- 
tribution presents the main objective of our Letter. Only 
upon ensemble averaging, the averaged value of subveloc- 
ity coincides with one given by solving analytically the 
fractional Fokker-Planck equation (FFPE) 0, [l] . 

We pose next the question about the universality class 



2 



of subvelocity distributions in this setup and find a simple 
result which we confirm with numerical simulations. It is 
astounding that this universality class does not depend 
on imposing a periodic static potential V{x) in addition 
to an applied constant force F (what seems feasible ex- 
perimentally). It also does not depend on temperature 
given that a is temperature-independent. This universal- 
ity feature is similar in nature with an established uni- 
versal scaling relation 0, [1| between anomalous current 
and biased diffusion, originally suggested for simple (i.e., 
in absence of a periodic potential) biased subdiffusivc 
CTRW transport [U, 01 . It is derived below by use of an 
heuristic argumentation, i.e. by reduction to a coarse- 
grained CTRW. On the level of ensemble-averaged quan- 
tities, we therefore obtain a universal law for the relative 
fluctuations of (sub-)velocity, or fluctuations of anoma- 
lous current, which can be tested experimentally. 

Theory for the biased CTRW. We start out from a 
CTRW on a one-dimensional lattice with period length a 
and WTD ^ . The walk is biased and the nearest neigh- 
bors jumps (this assumption is not restrictive and can 
be relaxed) occur with force-dependent splitting proba- 
bilities (7+ (towards the right) and q~ (towards the left), 
with q'^ + = 1. After n steps, the mean displacement 
is Ln = (x) = na(q^ — q^). From now on, we measure 
distance L in the units of a*{F) = a{q~^ — q~). Time will 
be measured in units of r and the subvelocity in units 
of voiF) = r(l + a)a*/r". 

In the fixed time setup, we fix the final time t and 
ask about the probability p{n, t) to make n steps. The 
answer is well-known, see p. 248 in Ref. ^2]: In the 
Laplace-domain, it reads 



p{n,u) = L^[^(u)]». 



(3) 



For w — > (i.e. for t — > oo), the leading term expansion 
of the Laplace-transform of WTD in Eq. ([1]) is ip{u) ^ 
1 - u". This leads to 

pin,u) ~u"-iexp[nln(l-u")] ~ u"-^ exp(-7iw") (4) 

in the limit of large n — > oo. The expression exp(— rtw") is 
related to the characteristic function of the extreme Levy 
stable law Ca(t) of index a scaled with parameter n, i.e. 
in the time domain it corresponds to n~^^"jCa 
Considering n as a continuous parameter (distance L = n 
in units of a*) and noting that Eq.Q equals the Laplace 
transform of 



* d j_ ( _t_ 



dL Li/c 



where Ca{x) is the cumulative distribution function of 
the extreme Levy stable law, i.e.. 



It ft 



Thus we can extract the distribution for the FT- 
subvelocity Va = r(l + a)L/t°' via a change of random 
variable from L to Va, yielding in terms of the scaled 
subvelocity Cq = Va/vo{F), for all F > 0, the stationary 
one-sided Levy-stable distribution: 
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(5) 



This universal form of the subvelocity distribution 
presents a major result of our study. 

Let us demonstrate that the very same result is re- 
covered also within the time- of- flight setup. That is, we 
are looking for the asymptotic distribution of times to 
make a large number of n steps. The corresponding dis- 
tance will assume a sharply peaked distribution around 
its mean value which can be identified with the sample 
size L. The random time t necessary to traverse the 
system of length L in the TOF setup is essentially the 
time necessary to make n steps. The overall time to 
make n 3> 1 steps tends in distribution to a one-sided 
Levy law n~^/°'jCa {n~^^"t). To see this it is sufficient to 
notice that the Laplace transform of the probability to 
find this time is given by p{u) — [V'(m)]" ~ (1 — u")" ~ 
exp[rtln(l — u")] — exp(— rtw"). The distribution of Va is 
then be obtained by the same change of variable as used 
above (the only difference being that t is now a random 
variable and L is fixed) to yield the same result in ([5]). 

The averaged value of the scaled subvelocity with dis- 
tribution (O is one, = 1, i.e. the subvelocity in ^ is 
scaled in fact through its averaged value Va {F) = vq (F) . 
All the higher moments can be obtained using the change 
of variable y — [r(l -I- a)/va]^^°' and the relation 



y'^Caiy)dy 



r(l -r?/a) 
r(l - ,/) 



(6) 



which is valid for any rj G (— oo, a), see in Ref. [l3|. With 
77 = —2a, it yields the second moment 



Sa 



2r2(l-f a) 



r(l-H2a) 

The relative fluctuation of subvelocity Svq 



(7) 



'a ~ {vaY /va equals in fact the universal scaling 
relation between the averaged subdiffusion current and 
biased diffusion of Refs. i-e-, 



we obtain, by a change of variable of integration from t 
to ^ = t/ii/". 



d r'/^"" d ( t 



5x^{t) 



/2r2(l + a) 



6x{t) 



r(l-|-2a) 
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t — ^00 



{6x{t)) 



(8) 
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This result is not trivial: This is so because (...) is 
the average over the stationary subvelocity distribution 
p{va) = p{Ca = Va/va)/va, whilc (....) is the average 
over the time-dependent population probabilities pi (t) of 
the lattice sites. This constitutes our second main result; 
it shows that weak ergodicity breaking is at the root of 
this remarkable scaling relation ([5]) . Namely, it is respon- 
sible for the startling change of the law of subdiffusion 
from {6x^{t)) oc when F = to {Sx^{t)} oc at 
F ^ 0, i.e. subdiffusion turns over into superdiffusion 
for 0.5 < a < 1. We remark, however, that the scaling 
relation ([5]) between the current and the biased diffusion 
cannot be employed to deduce the main result in ([5]). 

In particular, for a — 0.5 Eq. ^ simplifies to one- 
sided Gaussian form (cf. Fig. [T|) 



p(Ci/2) = -exp (--Cm) 



and 5vi^2 — ^7r/2 — 1. For other values of a, a handy 
approximation to the subvelocity distribution can be ob- 
tained using the corresponding small- a; asymptotic be- 
havior of the Levy-stable distribution [H, [3]. It reads, 



p(Cq) ^ A{a){aCa] 



exp[-B{a)C, 



l/(l-a)l 



(10) 



where A{a) = [i/27r(l - a)r(l + a)]-^ and B{a) = 
(1 - a)a"/(i-")r(l -I- a)-i/(i-"). For a = 0.5 this 
approximation reproduces the exact result in ([9|). For 
0.5 < a < 1, it correctly predicts that the distribution 
function is non-monotonic, possessing a maximum, cf. 
Fig. [2] below, which becomes sharp for a —^ I. In this 
limit, the relative fluctuation vanishes, Sva — * 0, and the 
velocity distribution tends to the delta-function centered 
at Va- However, the correct value p{0) always remains 
finite for a < 1, implying that there are always particles 
which become immobilized. For a < 0.5, p{va) decays 
monotonically. Moreover, for small a —^ 0, the distribu- 
tion becomes nearly exponential, consistent with Sva — > 1 
in this limit, see Fig. [3] below. 

All our analytical findings are confirmed by the nu- 
merical simulations of the underlying CTRW in a biased 
cosine potential V{x) = Vq cos{2ttx/1), using the Mittag- 
Leffler distribution ipiT) ^^'^ the numerical algorithm de- 
tailed in [1, . It is surprising that all these results hold 
universally, i.e. these are independent of the details of 
periodic potential and the temperature. To elucidate this 
astounding fact, being numerically confirmed with Figs. 
[U [21 E] for a washboard potential (details are given be- 
low) , we make use of the reasoning put forward with Ref . 
0. 

Theory for washboard potentials. We dilate the lat- 
tice by introducing many more points with separation 
Ax — > 0. The residence time distributions on every point 
are chosen to be Mittag-Leffler distributions Q with dif- 
ferent time scaling parameters = l/vi. This distribu- 
tion belongs to the same class in ([T|). Each point i is char- 
acterized also by the right and left nearest neighbor jump 
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FIG. 1: (Color online) Numerical subvelocity distribution 
Voi) for a — 0.5 in both FT and TOF setups for a peri- 
odic potential V{x) — Vo cos{2'kx/1) and differing bias forces 
F. Va is scaled in units of Vo{F) — Va{F)/vcr, where Va is 
given by (lllf) and Vcr ~ Fcri^a/{kBT). The solid lines depict 
the theoretical result ©I p{va) ~ p(Ca = Va /v*-,) /v*-,. 




FIG. 2: (Color online) Numerical subvelocity distribution 
p{va) vs. the analytical approximation in Eq. (|lUp and the 
exact result in Eq. ([5)l for a = 0.8. The same scaling applies 
as in Fig. [T] 



probabilities qf = gf /(g^ + g^), and by the fractional 
forward and backward rates, gf^ — q^vf, respectively. 
These quantities follow from the potential U {x) as so that 
the Boltzmann relation g^^i/g^ = exp[l3{Ui-i — Ui)] is 
fulfilled. Here, U, = U{iAx), C/,±i/2 = UiiAx ± Ax/2), 
and i^i — {gf + g^Y^"' ■ The generalized master equation 
for such a CTRW is [EQ: 



=5+1 P^-l{t) + g- 1 P,+,{t) - {g+ + gr)P{t) . 



In the spatial continuous limit Ax — > 0, it yields the 
FFPE (12). In this way, we simulate the stochastic dy- 
namics associated with on a sufficiently dense grid 
with step Ax, using the Monte Carlo algorithm from 0]. 



4 



10' 



10" 



_8 



10" 



10- 



1 1 


1 

FT o : 




TOF ^ 




Eq.(5) -■ 

Eq.(lO) ; 




F=l 




- 


V, V 


^ aOW O a 

A 

|o\ o 



FIG. 3: (Color online) Same as in Fig. [T]for q = 0.2. 



Consider next a periodic potential with period I sub- 
jected to a finite bias force F. One can course-grain 
the corresponding limiting CTRW and to map it onto 
a new biased CTRW with the lattice period I, i.e. we 
average over spatial period I. The precise form of the 
coarse-grained WTD is not known. However, it belongs 
to the same class as ([Ij ; only the time parameter r is cor- 
respondingly changed together with the coarse-grained 
splitting probabilities q^. We note that such course- 
graining of Markovian, normal diffusion in washboard 
potentials yields a non-Markovian CTRW which can give 
rise to such profound effects as giant acceleration of dif- 
fusion In clear contrast, our original CTRW is al- 
ready a profoundly non-Markovian, non-ergodic process 
possessing infinite memory. Coarse-graining it further 
does not change the universality class because no corre- 
lations between the residence times in non-overlapping 
spatial domains occur. 



For arbitrary periodic tilted potentials, the result for 
ensemble-averaged subvelocity was obtained in Refs. 0, 
i as: 



kJ [1 - exp{-f3Fl)] 



/Jdx/;+'dy exp(-/3[C/(x)-f/(2/)]) 



(11) 



In all our numerical simulations we used the archetype 
cosine potential V{x) — Vq cos{2nx/l). The grid contains 
200 points within each spatial period. A scaled temper- 
ature of ksT — 0.1 Vb is used throughout and the force 
F is scaled in units of the critical force F^r, where with 
F > Fcr the potential U{x) becomes monotonic without 
barriers in between. The number of particles is A'^ = 10^. 
The different lines for fixed a and different bias values 
F, are due to the different values of the scaling parame- 
ter Va{F), calculated in accordance with In accor- 
dance with our theory, all the related lines perfectly co- 
incide (not shown) after re-scaling Va Ca — Va/va{F), 
p ^ p ■ Va{F), for all F > 0. The numerical results thus 
corroborate with theory. 

In conclusion, we have shown that the weak ergodicity 
breaking is responsible for the universal scaling relation 
(III) between the anomalous current and subdiffusion oc- 
curring in arbitrary tilted periodic potentials. This in- 
triguing result follows from the universal law for the the- 
oretically deduced subvelocity distribution in ^ which 
is the major finding of this work. 
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